We report preliminary results of a search for CPT and Lorentz violation in 
B^-B^ oscillations using an inclusive dilepton sample collected by the BABAR 
experiment at the PEP-II B Factory. Using a sample of 232 million BB pairs, 
we search for time-dependent variations in the complex CPT parameter z = zo + 
Zi cos {Vlt + (/)) where $7 is the Earth's sidereal frequency and t is sidereal time. 
We measure Imzo = (-14.1±7.3(stat.)±2.4(syst.)) x IQ-^, AFxRezo = (-7.2± 
4.1(stat.)±2.1(syst.)) x lO-^pg-i, Imzi = (-24.0 ± 10.7(stat.) ± 5.9(syst.)) x 
10-3, and Ar X Rezi = (-18.8 ± 5.5(stat.) ± 4.0(syst.)) x lO-^ps"!, where 
Ar is the difference between the decay rates of the neutral B mass eigenstates. 
The statistical correlation between the measurements of Im Zq and AF x Re Zg is 
76%; between Imzi and AF x Re Zi it is 79%. These results are used to evaluate 
expressions involving coefficients for Lorentz and CPT violation in the general 
Lorentz- violating standard-model extension. In a complementary approach, we 
examine the spectral power of periodic variations in z over a wide range of 
frequencies and find no significant signal. 
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1 INTRODUCTION 



It has been shown [1] that "If CPT invariance is violated in an interacting quantum field theory, 
then that theory also violates Lorentz invariance." The general Lorentz-violating standard-model 
extension (SME) [2] has been used to show that the parameter for CPT violation in neutral meson 
oscillations depends on the 4- velocity of the meson [3]. In studies of T{4S) BB decays at 
asymmetric-energy e^e~ colliders, any observed CPT asymmetry should vary with sidereal time 
as the T{AS) boost direction rotates together with the Earth [4], completing one revolution with 
respect to the Universe in one sidereal day 0.99727 solar day). We report a search for such 
effects using inclusive dilepton events recorded by the BABAR detector at the PEP-II collider. 

The physical states of the B^-B^ system are eigenstates of a complex 2x2 effective Hamiltonian 
and may be written as 

\Bl) = p^/l^\B'') + q^/l + ^\B''), 

\Bh) = v^^^\B'')-q^^l^\&), (1) 

where L and H indicate "light" and "heavy." The complex parameter z vanishes if CPT is preserved. 
T invariance implies \q/p\ = 1, and CP invariance requires \q/p\ = 1 and z = 0. 

The leading-order CPT-violating contributions in the SME imply z depends on the meson 4- 
velocity P'^ = 7(1,/?) in the observer frame as [5] 



Am-iAT/2' 

Here Ao^ is the real part of the difference between the diagonal elements of the effective Hamilto- 
nian, and the magnitude of the decay rate difference AT = Th — is known to be small compared 
to the B^-B^ oscillation frequency Am = mn — fnL- The sidereal time dependence of z arises 
from the rotation of P relative to the constant vector Aa. The Aa^ contain flavor-dependent CPT- 
and Lorentz-violating coupling coefficients for the valence quarks in the B^ meson. Analogous, but 
distinct, Aa^ apply to oscillations of other neutral mesons. Limits on the Aa^ specific to K^K^ 
oscillations [6] and on the Aa^ specific to D^D^ oscillations [7] have been reported by the KTeV and 
FOCUS collaborations, respectively. KTeV has also reported constraints on sidereal-time variation 

of the CPT violation parameter 0_| [8]. 

We approximate the 4-velocity of each B meson by the T(4S') 4-vclocity so that z is common 
to each i? in a pair. We choose the meson 3-velocity to lie along —z in the rotating laboratory 
frame shown in Fig. 1. The non-rotating frame containing the constant vector Aa has Z along the 
Earth's rotation axis, corresponding to declination 90° in celestial equatorial coordinates. X and 
y, each in the equatorial plane, lie at right ascension 0° and 90°, respectively. At sidereal time 
t = z lies in the X-Z plane and y is coincident with Y . The CPT parameter z may then be 
expressed as 

z = z(i) = — . ,^ Aao — /3Aaz cos y — /3 sin v f Aay sin ili + Aax cos fit] , (3) 

Am — zAT/z L \ /J 

where cosx = z ■ Z and f2 = 27r/24rad-sidereal-hour^^ is the Earth's sidereal frequency. 

We use the latitude (37.4° N) and longitude (122.2° W) of the BABAR detector, together with 
the Lorentz boost of the T{AS) {(3^ = 0.55 directed 37.8° east of south), to determine cosx = 0-63 
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Figure 1: Basis {x,y,z) for the rotating laboratory frame, and basis {X,Y,Z) for the fixed non- 
rotating frame. The laboratory frame precesses around the Earth's rotation axis Z at the sidereal 
frequency i?. The angle between Z and the direction z opposite to the T(4S) boost direction at 
PEP-II is X = 51°. 



and i = {ic + 7.3) sidereal- hours, where to is Greenwich Mean Sidereal Time (GMST) for each 
event; hence 

1.14Aao - 0.35Aa^ - 0.43 f Aay sin Qt + Aax cos f2t) 
^<*' = ^ Am-iAT/2 W 

We convert the "timestamp" that records when the event occurred to the Julian date (J) and 
calculate GMST as specified by the U.S. Naval Observatory [9]: 

to = mod(18.697374558 + 24.06570982441908 D, 24), (5) 

where D = J — 2451545.0 is the number of Julian days since 12^: 00 Universal Time on January 1, 
2000. 

The clock used to set the event timestamp has a rate governed by the PEP-II 59.5 MHz master 
oscillator and is resynchronized with U.S. time standards via Network Time Protocol at intervals of 
less than four months. During such intervals the event timestamps are conservatively estimated to 
accumulate absolute errors of less than 30 seconds per month. The sidereal phase of each timestamp 
is therefore determined to better than 0.2%. 

Since sidereal time gains 12 hours every six months relative to solar time, possible day/night 
variations in detector response tend to cancel over sidereal time for long data-taking periods. The 
data used in this analysis were accumulated over a period of more than four years. 

Inclusive dilepton events, where both B mesons decay semileptonically {b Xii^, with i = e 
or /u), comprise 4% of all T(4S') — > BB decays and provide a very large data sample for studies of 
CPT violation in mixing. In direct semileptonic neutral B decays, the flavor B^{B^) is tagged by 
the charge of the daughter lepton ^+(£~). 
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At the T{4S) resonance, neutral B mesons are produced in a coherent P-wave state. The B 
mesons remain in orthogonal flavor states until one decays, after which the flavor of the other B 
meson continues to evolve in time. Neglecting second order terms in z, the decay rates for the three 
semileptonic decay configurations {i'^l'^, £~i~, £'^£~) are given by 

A^++ oc e-'^l^*l|p/gp{cosh(ArAt/2)-cos(AmAt)}, 
N— oc e-^l^*l|g/pp {cosh(ArAt/2) - cos(AmAt)} , 

N+- oc e-^l^*l {cosh(ArAt/2) - 2 Re z sinh( ArAi/2) + cos(AmAt) + 2Imzsin(AmAt)} ,(6) 

where T is the average neutral B decay rate, and At is the difference between the proper decay 
times of the two B mesons. The sign of At has a physical meaning only for opposite-sign dileptons 
and is given by At = t~^ — 1~, where t~^(t~) corresponds to £~^(£~), respectively. 

The opposite-sign dilepton CPT asymmetry Acpt, between events with At > and At < 0, 
compares the oscillation probabilities P{B^ B^) and P{B^ B^) and is sensitive to CPT 
violation through the parameter z: 

P(gO ^ ^0) _ p(^o ^ ^0) _ N+-{At>0)-N+-{At<0) 
^'^ ~ P(S0 ^ BO) + P(SO ^ 5°) ~ iV+-(At >0)+iV+-(Ai<0) 
2-Rezsinh(ArAt/2) + Imzsin(AmAt) 
~ cosh(ArAt/2) + cos(AmAt) ' ^ ' 

The experimental bound on |Ar| [10] is sufficiently small for the approximation Rezsinh(ArAt/2) ~ 
Ar X Rez X (At/2) to be valid over the range —15 < At < 15 ps used in this analysis, and we 
measure the product AT x Re z instead of Re z alone. 

Wc present measurements of Im z and AT x Re z using a simultaneous two-dimensional likelihood 
fit to the observed At and sidereal time (t ) distributions of opposite-sign and same-sign dilepton 
events. Inclusion of the same-sign events allows a better determination of the fraction of non-signal 
events (called "o6c" in Sect. 3) in which the lepton from one B meson is not a direct daughter. We 
search for variations in z of the form 

z = zo -I- zi COS {f2t + 4>) (8) 

with a period of one sidereal day, and extract values for the CPT- and Lorentz-violating coupling 
coefficients Aa^ in the SME from the measured quantities Imzo, Imzi, AT x Rezo, and Ar x Rezi. 
This extends our previous sidereal-time- independent analysis that measured Imzo, AT x Rezo, and 
\q/p\ with the same events [11]. In the decay rates, we use |Ar| = 6x 10^^ ps~^ in the cosh(ArAt/2) 
term and \q/p\ = 1, consistent with the values reported in Ref. [10] and Ref. [11], respectively. In 
a complementary approach, we use the periodogram method [12], developed for studies of variable 
stars, to detect directly any periodic variations in z over a wide range of frequencies and to measure 
their spectral power P(z^). 



2 THE BABAR DETECTOR AND DATASET 

This analysis is based on about 232 million T{AS) — > BB decays collected during 1999-2004 with 
the BABAR detector at the PEP-II asymmetric-energy e^e~ storage ring. An additional 16 fb^^ of 
"off-resonance" data recorded 40 MeV below the T(45') is used to model continuum background. 

The BABAR detector is described in detail elsewhere [13] . This analysis uses the tracking system 
composed of a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH), the 
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Cherenkov radiation detector (DIRC) for charged tt-K discrimination, the CsI(Tl) calorimeter 
(EMC) for electron identification, and the 18- layer flux return (IFR) located outside the 1.5-T 
solenoid coil and instrumented with resistive-plate chambers for muon identification and hadron 
rejection. A detailed Monte Carlo program based on GEANT4 [14] is used to simulate the response 
and performance of the BABAR detector. 

3 ANALYSIS METHOD AND LIKELIHOOD FIT 

The event selection is similar to that described in Rcf. [15]. Non-BS background, mainly due to 
e+e" qq {q = u, d, s, c) continuum events, is suppressed by applying requirements on the shape 
and the topology of the event. 

Lepton candidate tracks must have at least 12 hits in the DCH, at least one z-coordinate 
measurement in the SVT, and momentum between 0.8 and 2.3 CeV/c in the T(4S') rest frame. 
Electrons are selected by requirements on the ratio of the energy deposited in the EMC to the 
momentum measured in the DCH. Muons are identified through the energy released in the EMC, as 
well as the strip multiplicity, track continuity, and penetration depth in the IFR. Lepton candidates 
arc rejected if their signal in the DIRC is consistent with that of a kaon or a proton. The electron 
and muon selection efficiencies are about 85% and 55%, with pion misidentification probabilities 
around 0.2% and 3%, respectively. 

Electrons from photon conversions are identified and rejected with a negligible loss of efficiency 
for signal events. Leptons from J/ip and ip{2S) decays are identified by pairing them with other 
oppositely-charged candidates of the same lepton species, selected with looser criteria. The event 
is rejected if the invariant mass of any such lepton pair satisfies 3.037 < m£+£- < 3.137 GeV/c^ or 
3.646 < m£+£- < 3.726 GeV/c^. Remaining events with at least two leptons are retained, and the 
two highest momentum leptons in the T{4S) rest frame are used as the dilepton candidates. 

Separation between direct leptons {''^ b (. " ) and background cascade leptons from the " ^ 
c ^ £ " decay chain is achieved with a neural network that combines five discriminating variables: 
the momenta of the two lepton candidates, the angle between the momentum directions of the 
two leptons, and the total visible energy and missing momentum in the event, all computed in the 
T{AS) rest frame. Of the original sample of 232 million BB pairs, 1.4 million pass the selection. 

In the inclusive approach used here, the z coordinate of the B decay point is approximated by 
the z coordinate of the lepton candidate's point of closest approach in the transverse plane to our 
best estimate of the [x, y) decay point of the T{4S). In the transverse plane, both the intersection 
point of the lepton tracks and the beam-spot position provide information about the T(4S') decay 
point. We combine this information in a x^"fit that optimizes our estimate of the T[AS) decay point 
in the transverse plane using the transverse distances to the two lepton tracks and the transverse 
distance to the beam-spot position. The proper time difference At between the two B meson decays 
is taken as At = Az/{Pj)c, where Az is the difference between the z coordinates of the B decay 
points, with the same sign convention as for At, and (/?7) = 0.55 is the nominal Lorentz boost. 
For same-sign dileptons, the sign of At is chosen randomly. 

A large control sample of e^e~ ^'^^~{'^) events, with true Az = 0, was used to check for 
any sidereal-time-dependent bias in the Az measurement that could mimic a signal for Lorentz 
violation. The measured amplitude for such a bias at the sidereal frequency is (0.015 ± 0.025) //m, 
consistent with no variation around the mean value {Az) = (0.030 it 0.018) /xm. The corresponding 
amplitude for a sidereal-time-dependent bias in At for BB events is (9 ± 15) x 10~^ ps. Similar 
amplitudes are found for possible day/night variations in the Az and At measurements. 
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We model the contributions to our sample from BB decays using five categories of events, i, 
each represented by a probability density function (PDF) in At and sidereal time t, denoted by 
prijc fpj^gjj. shapes are determined using the B^B^ (n) and B'^B~ (c) Monte Carlo simulation 
separately, with the approach described in Ref. [16]. 

The five categories of dilepton BB decays, with contributions estimated from Monte Carlo 
simulation, are the following. Pure signal events with two direct leptons {sig), comprising 81% of 
the selected BB events, give information about the CPT parameter z. "Opposite B cascade" {ohc) 
events, where the direct Icpton and the cascade Icpton come from different B decays, contribute 
about 9%. "Same B cascade" {she) events, in which the direct lepton and the cascade lepton stem 
from the same B decay, contribute around 4%. About 3% of the dilepton events originate from 
the decay chain h ^ t~ ^ i~ " (Idlr), which tags the B flavor correctly. The remaining BB 
events [other) consist mainly of one direct lepton and one lepton from the decay of a charmonium 
resonance from the other B decay. 

The signal event PDFs, V^^g, are the convolution of an oscillatory term containing the sidereal- 
time dependent CPT parameter (Eq. 6) for neutral B decays (or an exponential function for charged 
B decays) with a resolution function that is the sum of three Gaussians (core, tail, and outlier) with 
means fixed to zero [17]. The widths of the narrower core and tail Gaussians are free parameters 
in the fit to data; the width of the outlier Gaussian is fixed to 8ps. The fractions of all three 
Gaussians are determined by the fit: the tail and outlier fractions are free parameters, and the sum 
of the three fractions is constrained to unity. 

The obc event PDFs, "P^fJ^, are modeled by the convolution of (At, t )-dependent terms, similar 
in form to those for signal, with a resolution function that takes into account the effect of the 
charmed meson lifetimes. Since both short-lived and D'^ , and long-lived mesons are involved 
in cascade decays, the resolution function for the long-lived and short-lived components is the 
convolution of a doublc-sidcd exponential with the sum of three Gaussians. To allow for possible 
outliers not present in the Monte Carlo simulation, the fraction of the outlier Gaussian is a free 
parameter in the fit to data. The parameterization of the she event PDFs, ■P"^^; account for the 
lifetimes of charmed mesons in a similar way. 

The PDFs for Idlr events, V^^^^^ are similar to those for the signal events. The resolution 
function takes into account the r lifetime and is chosen to be the convolution of two double- 
sided exponentials with two Gaussians. The PDFs for the remaining BB events, 'P^fher'' 
convolution of an exponential function with an effective lifetime and two Gaussians. 

The fractions (Z^^'^, /"ji,- and /"^'^g^) of she, Idlr and other events are determined directly from 

the B^B^ and B^B~ Monte Carlo simulations. The fraction /_| of B~^B~ events and the fraction 

of ohc events are free parameters in the fit to data. The ratio f obc/ f obc constrained to 

the estimate obtained from Monte Carlo samples. 

Non-i?i? events are estimated, using off-resonance data, to comprise fcont = (3.1 ± 0.1)% of 
the dilepton candidates. The PDF for this component is modeled using off-resonance dilepton 
candidates selected with looser criteria and on-resonance events that fail the continuum-rejection 
criteria. 

The CPT violation parameter z is extracted from a binned maximum likelihood fit to the events 
that pass the dilepton selection. The likelihood C contains the (At, t )-dependent PDFs described 
previously and combines 24 sidereal-time bins. 

The likelihood for each sidereal-time bin is given by 

'C(At) = fcontPcont + {i- fcont){f+-VB+B- f+-)'PsO^}^ "^^^^^ 

-p „-„ = (1 — f". )V'" ■+- f". "P" 
B^B^ \ J sig/'^ case' J sig'^ stg^ 



12 



' case — J other' other ' Jldlr' Idlr ' Jsbe' she ' Jobe' obc' 

Here we omit small charge asymmetries (~ 10^'^) induced by charge-dependent lepton reconstruc- 
tion and identification efficiencies. While affecting the same-sign decay rates, these asymmetries 
cancel at first order for opposite-sign dilepton events. Our previous sidereal-time-independent 
analysis [11] found these asymmetries to be a source of systematic uncertainty only for \q/p\- 

The likelihood fit gives Imzo = (-14.1 ± 7.3) x 10"^, AT x Rezo = (-7.2 ± 4.1) x 10"^ ps'^ 
Imzi = (-24.0±10.7)xl0-3, and AFxRezi = (-18.8±5.5) x 10"^ ps'^ The statistical correlation 
between the measurements of Imzo and AF x Rczq is 76%; between the measurements of Imzi 

and Ar x Rezi it is 79%. The fitted fractions of B^B^ and ohc events are /-) = (59.1 it 0.3)% 

and /^^^ = (10.7 ± 0.1)%, respectively. 

Figure 2 shows the asymmetry Aqpti defined in Eq. 7, as a function of sidereal time. The curve 
is a projection of the two-dimensional asymmetry onto the sidereal-time axis. To exhibit better the 
measured asymmetry, the projection is performed by integrating over |At| > 3ps thereby omitting 
highly-populated bins near |At| = where any asymmetry is predicted to be small and is diluted 
by At resolution effects. 
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Figure 2: The asymmetry Aqpt, integrated over |At| > 3ps, as a function of Greenwich Mean 
Sidereal Time in seconds folded over a period of 24 sidereal hours. The curve is a projection from 
the two-dimensional likelihood fit onto the sidereal time axis. 



4 SYSTEMATIC STUDIES 

There are several sources of systematic uncertainty in these measurements. To determine their 
magnitude, we vary each source of systematic effect by its known or estimated uncertainty, and 
take the resulting deviation in each of the measured parameters as its error. 

The widths of the core and tail Gaussians of the resolution function for the obc and sbc categories 
as well as the pseudo- lifetime for the Idlr and other categories are varied separately by 10%. This 
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variation is motivated by comparing the fitted parameters of the signal resolution function obtained 
from BB Monte Carlo samples and from data. The fractions of the short-lived and long-lived 
charmed meson components for obc and sbc are varied by 10%. Modeling of the PDFs is the main 
source of systematic uncertainty in ImzQ. 

We have also varied the parameters Am, r^o and r^i independently within their known un- 
certainties [18]. For Ar, we have allowed for 3a deviations from the value reported in Ref. [10] by 
varying |Ar| over the range 0-0.1 ps~^. The lifetimes r^o and r^i are the dominant sources of 
systematic uncertainty in Im zi and AF x Re zi. The dominant systematic uncertainty in AF x Rc zq 
is imperfect knowledge of the absolute z scale of the detector and residual uncertainties in the SVT 
local alignment. A possible sidereal-time-dependent bias in the At measurement with amplitude 
24 x 10~^ ps, derived from the amplitude (9 it 15) x 10~^ ps found with e'^e~ — n'^ix~{^) events, 
contributes a negligible systematic uncertainty. 

Table 1: Summary of systematic uncertainties for Imzo, AF x Rezo, Imzi, and AF x Rezi. 



Systematic Effects 


cr (Imzo) 


C7(AF X Rezo) 


(7(Imzi) 


cr(AF X Rezi) 




(xlO-3) 


(xlO-3 ps-i) 


(xlO-3) 


(xl0-3ps-^) 


PDF modeling 


±2.0 


±1.0 


±2.5 


±1.2 


Bkgd component fractions 


±0.1 


±0.1 


±0.2 


±0.2 


AF, Am, Tso, tb± 


±1.3 


±1.0 


±4.9 


±3.6 


SVT alignment 


±0.6 


±1.5 


±2.0 


±1.1 


Total 


±2.4 


±2.1 


±5.9 


±4.0 



For each parameter, the total systematic uncertainty is the sum in quadrature of the estimated 
systematic uncertainties from each source, as summarized in Table 1. 

5 CPT VIOLATION PARAMETER FREQUENCY ANALYSIS 

To perform a more general search for periodic variations in the CPT violation parameter z over 
a wide frequency range, we adopt the periodogram method [12] used in astronomy to study the 
periodicity of variable stars such as Cepheids. For each test frequency the method determines 
the spectral power P(i'), defined by 

1 , ^ 2 
n^)^J^\L^3^^''"^' ' (10) 



w 



J= 



from N data points measured at times Tj and having weights Wj with variance cj^. In our case Tj 
is the Universal Time of event j. In the absence of an oscillatory signal, the probability that the 
largest P(z^) exceeds a value S is given by 



Pv{P^^{u)>S;M} = l-(l-e-^) , (11) 



where M is the number of independent frequencies tested. 

Our search uses 20994 test frequencies from 0.26year~^ to 2.1day~^ in units of (solar time)~^, 
with steps of 10~^ day~^. To guard against underestimating the spectral power of a signal, we have 
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Figure 3: Periodograms for opposite-sign dilcptons showing spectral power P(i^) for weights wj oc 
Atj (top), Wj oc sin(AmAij) (center), and wj oc AmAtj — sm{AmAtj) (bottom) providing sensi- 
tivity to Ar X Rezi, to Imzi, and to Imzi subject to the SME constraint AF x Rez = 2AmImz, 
respectively. 
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oversampled the frequency range by a factor of about 2.2. The number of independent frequencies 
is about 9500. Twenty-seven test frequencies he between the Earth's sidereal and solar rotation 
frequencies. Each weight wj depends on the decay time difference Atj reconstructed for event 
j occuring at time Tj. Periodic variations in z affect the decay rate through the terms 

Imzsin(AmAt) and Rez sinh(ArAi/2) ~ AF x Rcz(At/2) in Eq. 6. Sensitivity to variations in 
Ar X Rez and Imz is attained by employing weights wj oc Atj and wj oc sin(AmAtj), respectively. 
In the context of the SME, the imaginary part of Eq. 2 implies AT x Rez = 2AmImz, and hence in 
Eq. 6 we have Imz sin (AmAt) — Rcz sinh (ArAf/2) ~ Imz [sin (Am At) — Am At ]. Accordingly, 
we also search for periodic variations in Imz using weights Wj oc AmAtj — sin(AmAtj). 

Figure 3 shows the spectral powers P(i^) measured in the opposite-sign dilepton data sample 
using the weights Atj, sin(AmAtj), and AmAtj — sin(AmAtj). The largest spectral power ob- 
tained for each of these weights corresponds to statistical fluctuation probabilities of 62%, 36%, 
and 76%, respectively, consistent with no periodic variation in the CPT violation parameter over 
the frequency range 0.26 year"-*^ to 2.1 day"-*^. At the Earth's sidereal frequency (r^ 1.0027 day^^), 
P(i/) = 3.73, 0.71, and 6.24 for the three weight types. At the Earth's solar-day frequency, the 
corresponding P(z^) = 1.50, 0.97, and 1.47. 

To check the validity of these results, we performed several tests of the periodogram method 
using events from data and from Monte Carlo simulation. Test periodograms showed large spectral 
powers at expected frequencies for (i) a generic dilepton Monte Carlo sample assigned event times 
Tj with a 0.5 sidereal-day ~^ frequency modulation, and (ii) unweighted dilepton data events, which 
give sensitivity to variations in the overall event rate — generally higher during night and weekend 
shifts, corresponding to frequencies of lday~^ and Iweek^^. Test periodograms for same-sign 
dilepton data events, which are not sensitive to CPT violation, showed no significant spectral 
power. The largest P(i^) value, obtained with Atj weights, corresponds to a statistical fluctuation 
probability of 7%. Test periodograms for opposite-sign dilepton data events, with the sign of At 
randomized to remove any measurable CPT violation, also showed no significant spectral power. 
We used these periodograms to check whether the distribution of P(z^) values has a probability 
density oc exp{— /c • P(j^)} with k = I, consistent with Eq. 11. A fit to the P(^') values yields 
k = 1.006 ± 0.001 with = 68.2 for 53 degrees of freedom. 

6 RESULTS 

Figure 4 shows confidence level contours for the parameters Imzi and AF x Rezi including both 
statistical and systematic errors. A significance of 2.2cj is found for periodic variations in the CPT 
violation parameter z at the sidereal frequency, characteristic of Lorentz violation. 

In the framework of the SME, the quantities Imzo, Imzi, AF x Rezo, and AF x Rezi are 
related by Eq. 4 to the Ao^ containing CPT- and Lorentz-violating coupling coefficients. With 
|AF| ^ Am, and using the SME constraint AF x Rez = 2AmImz implied by Eq. 2, we obtain 

1.14Aao - 0.35Aaz (Am/AF)AF x Rezo = 2Am(Am/AF)Imzo, 

(12) 

0.432 [(Aax)^ + (Aay)2] ^ [(Am/AF)AF x Rezi]^ = 4 [Am(Am/AF)Imzi]^ 

Taking into account error correlations between AF x Rezo and Imzo, and between AF x Rezi and 
Imzi, we find 

Aao - 0.30Aaz « -(5.2 ± 4.0)(Am/AF) x 10"^^ GeV, 
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Figure 4: Confidence level contours for the parameters Imzi and AT x Rezi including both statis- 
tical and systematic errors. The correlation between the measurements of Im zi and AF x Re zi is 
79%. The line contours indicate la, 2a, and 3a significance. The star at Imzi = AF x Rezi = 
indicates the condition for no sidereal-time dependence in z. 



(Aax)2 + (Aay)2 « (37 ± 16)| Am/AF| x 10^^^ GeV. 

Here we use Am = (0.507 ± 0.004) ps~^ = (3.34 ± 0.03) x 10"^^ GeV [19], and note that lattice 
QCD calculations give Am/ AT ~ —200 in the standard model [20]. 

7 CONCLUSIONS 

We have used data containing 232 million BB pairs to perform a simultaneous likelihood fit of 
same-sign and opposite-sign dilepton events that includes both the decay time difference At and 
the sidereal time t of each event. We have measured the CPT violation parameter of form z = 
zq + zi COS {f2t + (p) and find 

Imzo = (-14.1 ± 7.3(stat.) ± 2.4(syst.)) x 10"^ 

AFxRezo = (-7.2 ± 4.1(stat.) ± 2.1(syst.)) x 10"^ ps~\ 

Imzi = (-24.0 ± 10.7(stat.) ± 5.9(syst.)) x 10~^ 

AFxRezi = (-18.8±5.5(stat.) ±4.0(syst.)) X lO^^ps"^ 

A significance of 2.2a, compatible with no sidereal-time dependence, is found for periodic variations 
in z at the sidereal frequency that are characteristic of Lorentz violation. The complementary peri- 
odogram method provides no strong evidence for Lorentz and CPT violation, or for any periodicity 
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in z over the frequency range 0.26 year" to 2.1 day" . The results of the hkchhood fit are used 
to constrain the quantities Aa^^ containing CPT- and Lorentz-violating couphng coefficients for 
neutral B oscillations in the general Lorentz-violating standard-model extension. 
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